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ABSTRACT 

We have analyzed new HST/ACS and HST/WFC3 imaging in E475W and E814W of two previously- 
unobserved fields along the M31 minor axis to confirm our previous constraints on the shape of M3I’s 
inner stellar halo. Both of these new datasets reach a depth of at least F8I4W<27 and clearly detect 
the blue horizontal branch (BHB) of the field as a distinct feature of the color-magnitude diagram. We 
measure the density of BHB stars and the ratio of BHB to red giant branch stars in each field using 
identical techniques to our previous work. We find excellent agreement with our previous measurement 
of a power-law for the 2-D projected surface density with an index of 2.6j'jg'2 outside of 3 kpc, which 
flattens to a <1.2 inside of 3 kpc. Our findings confirm our previous suggestion that the field BHB stars 
in M31 are part of the halo population. However, the total halo profile is now known to differ from this 
BHB profile, which suggests that we have isolated the metal-poor component. This component appears 
to have an unbroken power-law profile from 3-150 kpc but accounts for only about half of the total 
halo stellar mass. Discrepancies between the BHB density profile and other measurements of the inner 
halo are therefore likely due to the different profile of the metal-rich halo component, which is not only 
steeper than the profile of the metal-poor component, but also has a larger core radius. These profile 
differences also help to explain the large ratio of BHB/RGB stars in our observations. 

Subject headings: galaxies: individual (M31) — galaxies: stellar populations — galaxies: evolution 


1. INTRODUCTION 


Stellar halos of massive galaxies are likely formed 
in great part from the disruption of dwarf galaxies. 
The stellar content of these halos thus o ffer quantita¬ 


tive in sight into galaxy a ss embly ( e.g. 


Johnston et al.l 

201(11 


200 


Deason et al.l 


Zolotov et al 

20131 


20011 Bullock fc JohnstonI 20051 Bell et all 


Pillepich et al 


Bullock et al.l 

20M 


2009 ; ICooper et al.l 


Im 


One 

key observable property that is typically compared to 
simulations is the radial stellar density profile. These 
profiles are comm only measured in nearby galaxies , 
such as the MW (iJuric et al.l l2008l ISesar et al. 1201111 


M31 (iGuhathakurta et ^ 2005 ~ ICourteau et al. 20111 


Ibata et ^ 2014b N253 (iBailin et ^20111 Greggio et all 

20141) and M81 ijBarker et al. 20091) . However, such pro¬ 


files have proven difficult to measure over wide radial 
ranges. In particular, the profile shape of the inner parts 
( ;$, 10 kpc) of stellar halos are difficult to probe. This 
shape is important to determine as the density profile at 
large radii can diverge if extrapolated to zero radius. 

Much work has gone into measuring the stellar density 
profile of M31. M31 makes an excellent case study because 
it is very nearby, and its halo has been resolved to 150 kpc 


fe.g. 

llbata et alJ 20011 Kalirai et alJ 2006al: Gilbert et al. 

2007 

: McConnachie et al. 20091 Gilbert et al.l 2009: 

Gourteau et al. 20111: Gilbert et al.l 2012: Williams et al. 

2012 

;lDorman et al.ll2013l: llbata et al.ll20l4 Gilbert et al. 

2014 

). Recently, the Panchromatic Hubble Andromeda 


Treasury (PHAT ibalcanton et al.l[2012ll has started a new 
era in the study of M31 by producing homogeneous HST 
imaging of about 1/3 of the star forming disk. One of 
the unexpected discoveries from this project was the de¬ 
tection of field BHB stars within M31 bulge fields. The 
surface density of these BHB stars was found to follow the 
trend from the outer halo, strongly suggesting that these 
stars belong to the halo population and allow us to directly 
measure the inner halo stella r density profile in M 31 to 
previously-unattainable radii (IWilliams et al.ll2012ll . 

Our work with the PHAT data suggested that, if the 
BHB is indeed a pure halo population, then the stellar 
halo has a surface density profile with a 2-D power-law 
index of a = 2.6jjn'o that flat tens to a <1.2 in the in¬ 
ner 3 kpc (IWilliams et al.l[2012fl , and the BHB/RGB ratio 
of the halo component is relatively constant with radius. 
These findings about the BHB population of M31 pro¬ 
duced predictions of the BHB stars that would be observed 
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in previously-unobserved HST fields. If confirmed, these 
would show that the power law continues to be steep to 
radii of ~3 kpc and turns over to > —1.2 at <3 kpc, the 
only measurement of a halo profile to very small radii. 

At the same time, there were a number of outstanding 
questions left open by the BHB stellar halo measurement. 
The interpretation of these BHB stars traced the halo com¬ 
ponent hinged on the assertion of a constant BHB/RGB 
ratio in the halo from 20 kpc (where it was measured in 
halo only fields) to the inner parts. Since the metallicity 
of the stellar halo is known to vary with radius, this asser¬ 
tion left open the possibility that inferences of the density 
profile of the inner stellar halo based on BHB stars were 
incomplete. 

Furthermore, recently there have been several rele¬ 
vant an d comprehensiv e stud ies of the M3 1 halo struc- 
ture by Dorm^an et al.l (l2013ll. Ibata et HI (|2014ll , and 
iGilbert et al.l ( 2014 1. iDorman et al.l ( 20131 1 decomposed 
the structure of M31 by simultaneously fitting surface 
brightness, stellar luminosity function, and stellar velocity 
data. They found a halo profile with a power-law index 
of 2.5±0.2, very close to that measured using the BHB 
stars alone. However, they found a core radius, where the 
halo stellar profile flattens, to be ~10 kpc, apparently at 
odds with the BHB stars, which clearly continue to fol¬ 
low the '^2.5 power-law index well inside of 6 kpc. Gon- 
firming previous measurements of a metallicity gradient 


in the M31 halo (|KaIiraLet al 


iTanaka et al.ll2010lb Ibata et al 


200^ iKoch et ai.l [20M 
( 2014) found that the dif¬ 


ferent metallicity components of the M31 halo have signif¬ 
icantly different profiles outside of 27 kpc, so that measur¬ 
ing the metal-poor component, as we trace with the BHB, 
does n ot necessarily tell the whole story. IGilbert et al.l 
(|2014ll quantified the metallicity gradient of the halo with 
spectroscopically-confirmed members and included radii as 
small as 9 kpc. 

In this paper, we probe BHB stars at intermediate radii, 
making it possible to test for departures from a ® 
power law suggested by the outer fields and the inner fields 
under the assumption that the BHB belongs to the halo 
component. We hnd that the BHB projected number den¬ 
sity matches a R“^ ® prohle to within the measurement 
uncertainties at these radii as well, strengthening the in¬ 
terpretation of the innermost BHB stars as probes of the 
metal-poor halo at small radii. This result combined with 
earlie r measurements of t he stellar halo d ensity at large 
radii (jGilbert et al.ll20l1 : llbata et al.ll2014D . yield a mea¬ 
surement of the metal-poor M31 stellar halo as an un¬ 
broken power-law with a 2-D index of ^2.6 over a very 
large radial range (3-150 kpc). Furthermore, the results 
suggest that the difference in the shapes of the metal- 
poor and metal-rich halo populations consists of differ¬ 
ent core radii as well as different power-law indexes. In 
section 2 we describe our new data and provide a brief 
reminder of the analysis technique. Section 3 discusses 
the results and comparison with predictions from previous 
work, and section 4 provides our conclusions. Through- 
out the paper, we assume a distance to M31 of 780 kpc 
(jStanek fc Garnavichl 1 199811. and a pos ition angle of the 
major axis of 38° ( Garignan et al.ir2006r) . 


2. DATA 


We observed the minor axis of M31 with HST for 
2 orbits, using the ACS as the primary camera and 
WFC3/UVIS as the parallel camera on 28 November 
2013. The telescope was oriented with PA_V3=259 de¬ 
grees, which allowed both the ACS and UVIS fields to fall 
on the minor axis simultaneously. The footprints of the 2 
fields are shown in Figure [T] Each camera obtained one 
orbit of data in F475W and one orbit in F814W. 

Image calibration was preformed by the MAST reduc¬ 
tion pipeline (OPUS version 2013_2) to produce the pro¬ 
cessed fit images on which we performed our photometry. 
The data were put through the same phot ometry pipeline 
as tha t used for the PHAT photometry in IWilliams et all 
(1201211 . making the reduction identical to that of our pre¬ 
vious study. Briefly, the fit images are put through mul¬ 
tidrizzle to identify cosmic rays; then they are multiplied 
by the pixel area map and bad pixels and cosmic rays are 
masked to prepare the images for photome try. Finally the 
images are put through the DOLPHOT (|Dolphinl l2000ll 
point spread function fitting routine to find and measure 
magnitudes for all of the stars in the image. DOLPHOT 
outputs quality metrics such as signal to noise sharpness 
and crowding, which are used to cull the results of artifacts 
and bad measurements. 

After this procedure 236561 and 61657 stars had been 
measured in the ACS and UVIS fields, respectively, with 
signal-to-noise >4, sharpness^<0.2, and crowding <2.25 in 
b oth bands. These qualit y metrics are discussed in detail 
in iDalcanton et al.l (l2012ll . The resulting color-magnitude 
diagrams are shown in Figure [21 where the smaller area of 
the UVIS camera as well as its slightly larger galactocen- 
tric distance cause the lower number of stars in the UVIS 
CMD. The BHB is clearly visible in both CMDs; however, 
the smaller number of BHB stars in the UVIS CMDs does 
result in somewhat larger uncertainties, especially in the 
fraction of BHB to RGB stars. 

Once the photometry was complete, we determined the 
number of BHB and R GB stars using the same method as 
IWilliams et al.l (120121) . In short, the number of BHB stars 
is measured by first fitting the F8I4W luminosity function 
in the color range 0.1<F475W-F814W<0.5 with the sum 
of a line and a Gaussian (see FigureE]). Once fitted, the in¬ 
tegral of the line alone (non-BHB stars) is subtracted from 
the integral of the fitted function to obtain the number of 
BHB stars. The number of RGB stars was measured by 
counting the number of stars in a box of CMD space where 
22.0<F814W<22.5 and 1.5<F475W-F8I4W<3.5. 


3. RESULTS AND DISCUSSION 


We have placed the measured BHB star densities and 
BHB/RGB fra ctions from ou r two new fields on the 
same plots as IWilliams et al.l (|20I2ll in Figures [4] and 
|5l The lines on the plots are not fits with the new 
data, but the same lines used in the previous analy¬ 
sis. The se lines correspond t o the halo profil e s mea ¬ 
sured by Williams et al.l (l20I2ll . ICourteau et al.l ( 2011 ). 
iGuhathakurta et ah ( 20051) . ^d Dorman et al.l ( 2013 ). 
The new obser vations are con s istent with the fits to 
the data from IWilliams et al.l (|2012D while probing a 
previously-unobserved radius. The fact that the density 
falloff mimics the expected drop from the inner radii to 
the outer radii shows that the BHB component of these 
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fields are not strongly affected by substructure. 

More importantly, the earlier data showed a steep tran¬ 
sition at these radii from low BHB/RGB ratio at the inner 
radii, where the RGB was highly contaminated by disk and 
bulge stars, to a relatively constant ratio of 0.7. The previ¬ 
ous measurements resulted in a prediction of a smooth in¬ 
crease in BHB/RGB ratio at these radii to the constant 0.7 
value. Indeed, the new observations confirm that the ra¬ 
tio follows the predicted smooth transition; however, new 
studies have since been published on the metallicity prop¬ 
erties of the halo with radius. In light of these studies, the 
interpretation of this transition is more complicated. 

Given that the BHB stars are tracing the metal-poor 
stellar halo, it is difficult to reconcile the fact that their 
density distribution clearly follows a continuous power-law 
to galactocentric distances inside of 6 kpc while sophisti¬ 
cated decompositions based on a much larger amount of 
data suggest that the stellar halo density distribution flat¬ 
tens outside of 8 kpc (|Gourteau et al.l[201ll : iDorman et al.l 
1201311 . 


Since our technique uses BHB stars, it necessarily only 
probes the metal-poor component of the stellar halo, sug¬ 
gesting that this component has a smaller core than the 
total stellar halo. Over the past decade, M 31’s halo has 
been shown to have a metallicity grad i ent (iKalirai et al.l 
2006bl: Koch et aP l2008t I Tanaka et al.l l201C)lb Recently. 


Ibata et al.l (1201411 found that the “unmaskecQ’ projected 


stellar density profile of the outer halo of M31 is a strong 
function of metallicity, with the more metal-rich compo¬ 
nent ([Fe/H]> —1.1) having a steep fall-off (a^3.7) at large 
radii. Their lower metallicity profiles ([Fe/H]< —1.1) have 
similar slopes to the BHB stars (2.3<a<2.7). Thus it is 
possible that the discrepancy between the BHB results and 
the kinematic and surface brightness results for the inner 
regions is attributable to the metal-rich and metal-poor 
components of the stellar halo having different profiles. 
This possibility affects the use of the BHB in understand¬ 
ing both the mass and profile of the M31 stellar halo. 

Below we explore the possibility that the explanation 
for the discrepancy between the multiple profile measure¬ 
ments is metallicity. Herein, we simplify the problem by 
breaking the halo into 2 components: one metal-rich and 
one metal-poor. M31’s halo is likely more complex than 
this simple two-component model— the distribution of 
metallicities in present in the population with radius is 
a mixture of stars with a range of metallicities. However, 
the BHB are likely to mainly represent the metal-poor 
populations in the h alo. As shown below , the large core 
radius measured by iDorman et al.l l|2013l) and the small 
core radius measured using the BHB stars can be some¬ 
what reconciled if (1) the observed differences between the 
profiles of the metal-rich and metal-poor stars persist into 
the inner regions of the halo, and (2) the large core ra¬ 
dius measured by Dorman et al. is comprised primarily of 
metal-rich halo stars. 


3.1. The Stellar Halo Profile 

Recent work confirms that the power-law profile of the 
metal-po or halo continues all the way out to at least 
150 kpc (|Gilbert et al.ll20l'^ llbata et alJl2014D . Thus, our 
mass estimate is valid for the metal-poor ([Fe/H]<-0.7) 

® measured without applying their masks to the substructure 


component of the halo. In particular, llbata et al.l (|2014li 
find that ^50% of the halo stellar mass from 27-150 kpc 
to have [Fe/Hl<— 0.6 (see their Table 4). The fraction is 
comparable to the lBrown et al.l[20081 SFH result. Clearly, 
a large amount of the halo stellar mass is metal-rich, which 
could change the characteristics of the core radius for the 
total halo. Thus, while the metal-poor halo appears to 
have a consistent power-law index of ^—2.5 from 3 to 150 
kpc, the total halo is more complex, with multiple power- 
law indices, a larger core, and a larger mass. 

If the large core radius of the total stellar halo is the 
result of the high-metallicity component, then one might 
expect the BHB/RGB ratio of the halo component to in¬ 
crease at smaller radii as the ratio of metal-poor stars to 
total stars increases. In FigurelHlwe plot the expected ra¬ 
tio of metal-poor stars to total stars with galactocentric 
distance. We assume that the metal-poor component of 
the halo has a profile with a core radius of 3.5 kpc and a 
2D power-law index of 2.7, and the metal-rich component 
has a core radius of 10.6 kpc and a 2D power-law index 
of 3.9. The result is that the metal-poor fraction has a 
minimum at 10-20 kpc, and a steep increase toward the 
galaxy center . Thes e ratios are similar to those found by 
iGilbert et al.l (l2014l) . which were ~0.2-0.7 between 9 and 
30 kpc, respectively. We are not able to match the spectro¬ 
scopic ratios exactly because the difference in power-law 
index between the two components does not result in the 
fraction changing that steeply. We note that the spectro¬ 
scopic sample in this ra dial range contains a significant 
amount of substructure (IGilbert et al.l 120071) . which may 
be causing the rapid change in the ratio of low-metallicity 
stars to total stars in those measurements. Therefore, we 
set the relative normalization between the two halo compo¬ 
nents so that in this radial range we roughly agree with the 
spectroscopic ratio of low metallicity stars to total stars. 

If we assume a BHB/RGB fraction of ^5, simila r to low- 
metallicity globular clusters (|Williams et al.l[2012fl , for the 
metal-poor component, and a fraction of 0 for the metal- 
rich component, then we can convert the metal-poor frac¬ 
tion into an expectation for the BHB/RGB ratio of the 
halo as a whole with galactocentric distance. This relation 
is shown in Figure [SI and roughly follows the metallicity 
fraction. 

Using t hese fractions along w ith the disk and bulge pro¬ 
files from IDorman et ^ (|2013f) . we calculate the expected 
BHB/RGB ratios for the total population, which are plot¬ 
ted in Figure |SJ This plot assumes that the metal-poor 
halo is the only component that has a significant BHB, 
and all components contribute to the RGB. The dashed, 
dotted, and gray lines assume a c onstant BHB/RGB hal o 
fraction of 0.7 and the pro files of Courteau et al.l (l201lD : 
iGuhathakurta et (l2005li . and IDorman et al.l (|2013ll . re- 
spectively. The solid line assumes the BHB/RGB fraction 
of the halo component follows the relation shown in Fig¬ 
ure |6l_and_the_totel profile follows that of the median fit 
from IPorman et "mI (l2013l) . Although still not a perfect 
match to the data, this relatively simple attempt to ac¬ 
count for the changing metallicity of the halo moves the 
model prediction considerably closer to the data in the in¬ 
ner regions. Therefore, one way to reconcile the measure¬ 
ments is through an increasing metal poor fraction inside 
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of the iDorman et al.l (|2013f) core radius. 

Therefore, one explanation for the difference between 
the core radii of the BHB profile and the total halo pro¬ 
file could be that the metal-poor halo has a significantly 
smaller core than the metal-rich halo. Such a possibil¬ 
ity may be testable in simulations, as many of the oldest, 
most metal-poor stars in the halo may actually form early 
and largely in the center of the dark matter halo. At the 
same time, if much of the metal-rich halo is formed later 
through accretion events, it could perhaps be hotter, re¬ 
sulting in a larger core. Our result can be compared to 
future simulations in this way. 

3.2. The Stellar Halo Mass 

If there is a significant amount of mass in the more 
metal-rich component of the halo which is not probed 
by the BHB stars, our previously-determined stellar halo 
mass of 2.ltJ;4XlO® Mq was not properly calibrated, 
as the profile was normalized to the total stellar sur- 
fa ce brightness at 21 kpc (^y'^29) from the photometry 
of iBrown et al.l (|2008ll . Based on observations of globu¬ 
lar clusters as well as stellar evolution theory, the BHB 
traces the metal-poor component of the stellar popula- 
tion. Accord i ng to the star formation history (SFH) in 
iBrown et al.l (l2008ll . about 50% of the stellar popula¬ 
tion has [Fe/H] -0.7. If we assume a H-band stellar 
mass-to-light ratio of of 2 for the metal- poor component 
(|Bruzual fc Char 1^120031 : IChabrieill2003ll . and normalize 
by the fraction of metal-poor stars (0.5), our metal-poor 
halo stellar mass measurement is 1.5^03X10® Mq. Since 
this only includes the metal-poor component, it is a lower- 
limit on the total stellar mass of the halo. 

We can move beyond a lower limit if we assume the 
metal-ric h component h a s a co re radius of 10 kpc, as mea¬ 
sured by iDorman et al.l (|2013ll and the ste e p 2-D power- 
law slope of -3.72 measured bv llbata et al.l (|2014ll for the 
unmasked metal-rich component. If we assume a H-band 
stella r mass-to-light ratio of 3 for the metal-ri ch compo¬ 
nent (|Bruzual &: Charl^ l2003t IChabrieJ l2003ll and nor¬ 
malize the metal-rich component to comprise the other 
half of the surface brightness at 21 kpc, then the mass 
of this component is ^1.4x10® Mq. Adding this to the 
mass of the metal-poor component then yields a total of 
^3x10® Mq. This estimate is similar to our earlier es¬ 
timate which did not properly account for the metal-rich 
component. Th is mass is a b out a factor of 2 lower than 
the estimates of llbata et al.l (|2014ll . Their mass estimates 
do not assume a mass-to-light ratio, but instead fit the ob¬ 
served stars to correct to the full underlying populations FI 
However, if we apply our assumed mass-to-light ratios of 


2 and 3 directly to their measured H-band luminosities 
for the metal-rich and metal-poor components of the outer 
halo, the masses of these components are 2 x 10® and 4x10® 
Mq, respectively. Summing these values and extrapolat¬ 
ing inward with an assumed core radius of 10 kpc yields 
an estimated total halo mass of ^2x10®. Thus, when the 
same assumptions are applied to both measurements, our 
halo stellar mass estimates are similar. However, ours is 
somewhat greater, largely due to the smaller core radius 
of the metal-poor component. 

4. CONCLUSIONS 

We have tested the predictions for the distribution of 
BHB stars from 30-40 arcmin (6-10 kpc) away from the 
center of M31 along the minor axis. The latest observa¬ 
tions confirm the predictions based on a model where the 
BHB stars trace the metal-poor M31 stellar halo, and the 
metal-poor stellar halo has a structure that follows a power 
law with 2-D index a = 2.6^q 2 inward to 3 kpc, then be¬ 
comes significantly flatter near the galaxy center. This 
core-radius is signihcantly smaller than that measured for 
the total stellar halo. The BHB/RGB ratio we measure 
inside of 10 kpc suggests that the metal-poor halo com¬ 
ponent may become increasingly importan t inside of the 
total h alo core radius of 10 kpc measured bv iDorman et al.l 
(|2013ll . Therefore, not only does the slope of the stellar 
halo profile change with metallicity, the core radius ap¬ 
pears to change with metallicity as well. 

Combining this result with recent measurements of the 
total halo and metal-ric h halo, suggests that the stellar 
halo mass calculated in (|Williams et al.l [201211 was likely 
underestimated because it did not properly account for the 
fact that the BHB stars were only probing the metal-poor 
component of the stellar halo. A more realistic calculation 
uses the metal-poor component, assuming a mass-to-light 
ratio of 2, to set a lower-limit on the stellar halo mass of 
1.51q 3 X 10 ® Mq. If we estimate the mas s of the metal-rich 
component using the llbata et aP (l2014ll profile, assuming 
a metal-rich mass-to-light ratio of 3, and normalizing to 
the surface brightness at 21 kpc, then we can add this 
to the metal-poor component estimate to derive a total 
stellar mass of the M31 halo is ^3x10® Mq. The signif¬ 
icant differences between the metal-rich and metal-poor 
halo characteristics in M31 should be reflected in galaxy 
formation and evolution simulations. 

Support for this work was provided by NASA through 
grant GO-12997 from the Space Telescope Science Insti¬ 
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for Research in Astronomy, Inc., for NASA, under contract 
NAS 5-26555. 
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Fig . 1.— Black rectangle s mark the location of our new ACS and UVIS fields along the minor axis of M31. Cya n rectangles ma r k the fields 
from [Williams et al.l (120121) that resulted in a BHB density measurement. While rectangles mark the fields from [Williams et ahl (120121) that 
allowed BHB density upper limit measurements. 
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Fig. 2.— Color-magnitude diagrams (CMDs) from the ACS {top) and UVIS {bottom) imaging data. The BHB is outlined. The UVIS field 
has significantly fewer stars because it is smaller and farther out. 
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Fig. 3.— Fits to the luminosity function near the BHB in our two new HST fields. A clear peak is detected, and the number of stars in the 
feature is given in the upper-left for each field. 
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Fig. 4.— The n umber density o f BHB stars meas ured for these fields (black data points at 30' and 36 ') overplotted on the radial 
distribution found in Williams et al.| ( gray points : 1201^. Solid, dashed, and dotted lines mark the halo fits from [Williams et al.l II2012D . from 
ICourtea u et al. l 112011 ), an d icThathakurta et al.l II2005I) . respectively. Light gray dot-dashed lines mark a random draw of 256 MCMC results 
from lPormai^t al.l l|2013l ). The new data in the previously-unobserved region agrees with the model profile determined assuming that the 
BHB traces the low metallicity component of the stellar halo. 
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Fig. 5.— The BHB/RGB ratio as a function of radius. Gray points are from IWilliams et al.l l|2012n . Black poi nts are from the new 
observations at 30' and 36'. Top: Dashed, and dotte d lines assume a consta nt BH B/RGB ratio of the pure h alo of 0.7 llWilliams et al.ll2012D 
and the halo fractions as a function of radius from ICourteau et al.l ll2011f) . and IGuhathakurta et ahl (120051), respectively. Gray lines also 
assume a constant BHB/RGB ratio of the pure halo of 0.7 but the halo fractions come from 256 trials from the lDorman et (120131) analysis. 
The strong transition between dominant components was predicted to occur in the previously-unobserved region from 20' to 60', and indeed, 
the new observations (at 30' and 36') follow the predicted sharp transition. Bottom: To match the sharp transition in BHB/RGB ratio, as 
well as the inner data points, requires an increase in the fraction of metal- poor halo stars in side of 10 kpc, as shown with the black solid 
line. This li ne is the result of extr apolating of the metal-rich halo profile of [lbata et al.l (120141) inward, adopting the stellar halo core radius 
measured bv IDorman et ahl (j2013D for the metal-rich component core radius, and adopting the BHB core radius and power-law index for the 
metal-poor halo component. Even this two-component halo model does not fully match the data; however, it is a significant improvement 
over any single-metallicity option, assuming previously-measured decomposition parameters. 
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F ig. 6. — Left: Th e fraction of metal-poor ([Fe/H]<-0. 7) halo stars a s a fu nction of galactocentric distance predicted by extrapolating 
the llbata et al.l ||2014I ~) metal-rich halo profile, adopting the lDorman et ^ HMH) core radius for the metal-rich component, and adopting the 
BHB profile for the metal-poor component. Because of the smaller core radius of the metal-poor component, the metal-poor fraction recovers 
inside of the metal-rich core radius. Right: The predicted BHB/RGB ratio of the halo alo ne as a function of gal actocentric distance, assuming 
the metal-poor component has a ratio of 5, comparable to metal-poor globular clusters ijWilliams et al.ll2Q12f) . The high value in the inner 
galaxy is again due to the smaller core radius of the metal-poor component. 












